Comparison between a diagrammatic theory for the BCS-BEC crossover and 

Quantum Monte Carlo results 
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Predictions for the chemical potential and the excitation gap recently obtained by our diagram- 
matic theory for the BCS-BEC crossover in the superfluid phase are compared with novel Quantum 
Monte Carlo results at zero temperature now available in the literature. A remarkable agreement 
is found between the results obtained by the two approaches. 
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The recent experimental realization of the BCS-BEC 
crossover with ultracold trapped Fermi atoms 1 has given 
impetus to theoretical investigations of this crossover. In 
a recent paper 2 , the t-matrix self-energy approach (orig- 
inally conceived for the normal phased) was extended to 
the superfluid phase, aiming at improving the description 
of the BCS-BEC crossover by including pairing fluctua- 
tions on top of the BCS mean-field approach considered 
m Refs.iandU 

In this theory, the effects of the collective Bogoliubov- 
Anderson mode is explicitly included in the fermionic 
self-energy, thus generalizing the theory due to Popov 
for a weakly- interacting (dilute) superfluid Fermi gas 6 . 
The theory is based on a judicious choice of the fermionic 
self-energy, such that it reproduces the fermionic mean- 
field BCS behavior plus pairing fluctuations in the weak- 
coupling limit as well as the Bogoliubov description for 
the composite bosons which form in the strong-coupling 
limit. In the intermediate-coupling region of interest 
about the unitarity limit, where no small parameter ex- 
ists to control the many-body approximations, the theory 
is able to capture the essential physics of the problem, as 
the excellent agreement with a previously available QMC 
calculation! at the unitarity point (kpap) 1 = has al- 
ready shown£, and as more extensively demonstrated by 
the present comparison with more recent QMC data-*i°" 
spanning the whole crossover region. The theory ofRef. 2 
is completely ab initio and it contains no adjustable pa- 
rameter. Although the comparison with QMC data is 
here limited to the zero-temperature limit where they 
are available, the predictions of the theory of Ref. ex- 
tend as well to finite temperature and across the critical 
temperature. 

Purpose of this Brief Report is to compare the the- 
oretical predictions obtained from the theory of Ref. 
with novel Quantum Monte Carlo (QMC) dat a 9 ' 10 , which 
were published after completion of Ref. |2J. A quantitative 
comparison between the results for the density profiles 
obtained from a local density version^ to the theory of 
Ref.0 and the experimental data was already presented 
in Ref. [H 

Both our calculations and the QMC calculations of 
Refs. @ and ^3 are based on a model Hamiltonian de- 
scribing a system of fermions mutually interacting via 
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FIG. I: Chemical potential at zero temperature vs the cou- 
pling parameter (feap)" 1 . The results of the present theory 
(t-matrix-I) and of its version without the inclusion of the 
self-energy shift Eo (t-matrix-II) are compared with the BCS 
mean field (BCS), the Fixed-node QMC data from Ref. [13 
(FNQMC), the Galitskii's expression for the dilute Fermi gas 
(Galitskii) , and the asymptotic expression for strong coupling 
using the result as = 0.6(jf. 



an attractive contact potential. In Ref. [T2| this Hamilto- 
nian was proved appropriate to describe the BCS-BEC 
crossover with trapped Fermi gases. In this model, the 
only dimensionless parameter representing the effective 
coupling strength is the (inverse of the) product kpap 
between the Fermi wave vector kp and the fermionic 
scattering length a p . For the homogeneous gas here con- 
sidered, kp — (37T 2 ?!) 1 / 3 where n is the particle density. 
Comparison will be made at zero temperature only, since 
finite-temperature QMC calculations for the BCS-BEC 
crossover are not yet available. 

The overall agreement between the two alternative (di- 
agrammatic and QMC) calculations turns out to be quite 
good, expecially in the most interesting intermediate- 
coupling regime about (kpap)^ 1 = 0. Figure 1 shows 
the comparison for the chemical potential at zero tem- 
perature, as obtained by our calculation^ and bv the 
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Fixed Node QMC (FNQMC) calculations of Ref. Efl. As 
discussed in Ref. |2j, on the weak-coupling side we find 
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FIG. 2: Excitation gap A m at zero temperature vs the cou- 
pling parameter (fepOF) -1 . The results of the present theory 
(t-matrix-I) are compared with the Green's function QMC 
data of Refs. and @ (GFQMC) as well as with the BCS 
mean field (BCS). 



it appropriate to introduce a constant shift So in the 
bare Green's function entering the self-energy. This shift 
needs to be included only for coupling values (kpap) -1 < 
—0.5, such that the self-energy can be considered to be 
approximatively constant. The curve obtained by this 
procedure is reported in Fig. 1 with the label t-matrix-I 
and corresponds to the data reported in Fig. 6 of Ref. 
For completeness, we also report in Fig. 1 the curve ob- 
tained without the inclusion of the self-energy shift So 
[with the label t-matrix-II] (by definition, the two curves 
I and II coincide when (fcp-Oi?) -1 > —0.5). Our results 
are in excellent agreement with the FMQMC data in the 
range —0.5 < (kpap) -1 ^5 0.5 spanning the crossover 
region. 

For couplings {k F a F y x < -0.5, the FNQMC results 
are extremely close to both our curves, lying just in be- 
tween them. In the weak-coupling region (fepcii?)" 1 < 
—2, our curves (as well as the FNQMC data) approach 
the asymptotic expression by Galitskiii^ for the chem- 
ical potential of a dilute Fermi gas. The BCS mean 
field (also reported in Fig. 1) misses instead the Gal- 
itskii correction to the non-interacting chemical poten- 
tial. More specifically, we have verified that our the- 
ory with the inclusion of the self-energy shift So [t- 
matrix-I] recovers the complete Galitskii's expression 
[ijtF — 1 + -^kpap + jjpr(ll — 21n2)(fc^ of) 2 including 
the second-order correction in kpap. The curve for the 
chemical potential obtained without the inclusion of the 
shift So recovers instead only the leading order correc- 
tion linear in kpap. [It can also be shown that neglect- 
ing the shift So introduces a spurious additional term 
^(-^kpap) 2 to the second-order correction in the Galit- 
skii's expression.] 

On the strong-coupling side, for coupling values 
(kpap) -1 > 0.5 our results deviate somewhat from the 



FNQMC data. This discrepancy is due to the fact that 
in our approach the boson-boson scattering is treated 
at the level of the Born approximation, corresponding 
to the value as = 2of of the bosonic scattering length 
a b ■ The importance of including the correct value of the 
bosonic scattering length (as = 0.6of, as calculated in 
Ref. IbH) in this region is clearly seen from the agreement 
between the FNQMC data and the asymptotic expression 
fx = — e /2 + /is/2, where eo is the binding energy of the 
2-body problem and [is = iTmsas / m B, with ng = n/2, 
m,B = 2to, and as — O.Gclf- The asymptotic curve cor- 
responding to the value as = lap almost coincides with 
our curve in this region. [This curve is not reported in 
Fig. 1 for overall clarity.] It is, finally, interesting to 
mention that the inclusion of the next-order correction 
to the bosonic chemical potential, corresponding to the 
expression fi B = 47r " B " B [1 + ^-Ubo^/it)] 1 ^ obtained 
in Ref. 15, would worsen appreciably the comparison be- 
tween the QMC data and the asymptotic curve in the 
coupling region 0.2 < (kpap)^ 1 ^, 2. The inclusion of 
this next-order term improves the comparison only in 
the truly asymptotic regime for {kpap) 1 ^ 2 (not re- 
ported in the figure) , where the next-order correction to 
the bosonic chemical potential is, however, already quite 
small. This finding could (at least partially) explain the 
absence of beyond-mean-field corrections on the bosonic 
side of the BCS-BEC crossover, recently reported in ex- 
periments with ultracold Fermi gasesi&. 

Quite generally, any theory of the BCS-BEC crossover 
connects the equation for the chemical potential jj, to the 
equation for the gap (order) parameter A in the super- 
fluid phase. The latter quantity is not directly accessible 
to the QMC simulations of Refs. and 0. In Ref. 0, 
however, the even-odd staggering of the ground-state en- 
ergy for a system with a finite number of particles was 
exploited to calculate the single-particle excitation gap 
A m . In a BCS-like framework (and for a sufficiently 
large number of particles) the gap A m is expected to 
coincide with the gap (order) parameter A when fi is 
positive and with the quantity (A 2 + /i 2 ) 1 ^ 2 when fj, is 
negative. For a given coupling, this gap occurs at the 
wave vector |k| = ^/2jl for positive fj, and at k = for 
negative /1. These BCS-like results are not expected to 
hold exactly away from weak coupling. The calculations 
presented in Ref. |2j, nevertheless, show that the identifi- 
cation of the single-particle excitation gap A m with A for 

> and with (A 2 + /1 2 ) 1 / 2 for [i < works fairly well 
for all couplings of interest. In particular, in Fig. 14 of 
Ref. this definition of the excitation gap A m was com- 
pared with the results obtained from an accurate analysis 
of the single-particle spectral function A(k, to), showing 
that the two definitions are in good agreement with each 
other over a wide coupling range. 

In Fig. 2 we compare A m , as obtained from our results 
for A and [i, with the QMC data of Ref. 9. The BCS 
mean field results are also reported for completeness. For 
the coupling value (kpap) -1 = a single QMC datum 
previously available from Ref. is also reported in the fig- 
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ure (full square). Even for the excitation gap, our results 
appear to be in remarkable agreement with QMC data 
in the crossover region —1 < (fci?ai?) _1 < 0.4. At larger 
couplings, the QMC results start instead to deviate from 
our results, the discrepancy being mainly due to the fi- 
nite range of the interaction potential used in the QMC 
calculations. In strong coupling, both our excitation gap 
and that calculated from QMC simulations tend, in fact, 
to half the value of the binding energy eo of the two-body 
problem. The binding energies for the contact potential 
and for the finite-range potential used in Ref. are close 
to each other only in a narrow range about (kpap) 1 = 0. 
At the coupling value (kpciF) 1 = 1, the binding energy 
for the finite-range potential of Ref. is already larger 
by about 40% than the contact-potential binding energy. 
This difference is responsible for the discrepancy between 
our values and the QMC data of Ref. on the strong- 
coupling side, where the excitation gap is controlled by 



the binding energy of the two-body problem. 

In conclusion, the theory of Ref. H for the BCS-BEC 
crossover in the broken-symmetry phase has been shown 
to compare extremely well with recent QMC data at 
zero temperature, especially in the intermediate-coupling 
(crossover) region which is the most interesting one both 
theoretically and experimentally. This agreement sug- 
gests that the choice of the fermionic self-energy made in 
Ref. captures the essential physics of the problem, as 
soon as the fermionic degrees of freedom get progressively 
quenched while forming composite bosons. 
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